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We present preliminary results for the decay constants of heavy-light mesons using NRQCD heavy and tadpole 
improved Clover light quarks. A comparison is made with data obtained using Wilson light quarks. We present 
an analysis of the 1/M corrections to the decay constants in the static limit and compare with the predictions of 
HQET. 



1. Introduction 

We present an analysis of the mass dependence 
of the heavy-light (HL) meson decay constants. 
We use NRQCD for the heavy quark and include 
all 1/M corrections to the static limit. The cor- 
responding NRQCD action is 
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In order to calculate matrix elements, and thus 
decay constants, to this order requires several 
operators. At tree level, these are obtained 
by relating the 4-component heavy quark field 
in full QCD (qh) to the 2-component NRQCD 
field (Q) using the Foldy-Wouthuysen transfor- 
mation (see [1]). Beyond tree level all opera- 
tors with the same quantum numbers can mix 
under renormalisation. Thus to 0(a/M) the ba- 
sis of operators for the axial-pseudoscalar current, 
975 To 9ft is given by ^ 
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Using translation invariance, currents involving 
O2 and O3 are identical on the lattice so it is only 
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necessary to calculate matrix elements of the tree 
level operators Oi and O'z- Similarly, the basis of 
operators for the vector current, qjiQu, 
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is reduced, for the numerical computation, to Oi, 
O2 and O3, where 7 matrix properties have also 
been used. 

In the same way as for the NRQCD action, the 
HL decay constants can be parameterised in the 
heavy quark limit in terms of an expansion in 
1/M. To 0(1/M) 
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Note that (7) is valid for both the PS and V decay 
constants; we define fVMy = {V \V \0) / ^/M^ . 
The coefficient cp is determined by nonper- 
turbative contributions arising from the hy- 
perfine interaction (Ghyp) and the kinetic en- 
ergy of the heavy quark (Gkin), which appear 
in Cnrqcd, and the corrections to the cur- 
rent (Gcorr cx (02 1 -PS') at tree level for /ps): 
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where rfjvf = 3 and —1 for PS and V mesons, 
respectively. From (8) the 0(1/M) coefficient 
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Similarly, the ratio of decay constants. 



is 1 in the 
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the 0(1/ M) coefficient Cp is given by 
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Thus, by computing various combinations of the 
PS and V currents the contributions from each 
0(1/M) term in the NRQCD action and matrix 
elements can be obtained. HQET provides an 
analogous decomposition of these coefficients (see 
[2] for a derivation of (8)). However, since HQET 
constructs an effective theory in terms of the 
heavy quark pole mass, as opposed to the bare 
quark mass in NRQCD, only physical combina- 
tions of Gi, ie Cp, Cp and Gkin, can be com- 
pared. Naively, these coefficients are expected to 
be O(Aqcd) ~ 200-500 MeV and the corrections 
to the static limit at Mb 0(Aqcd/Mb) ~ 10%. 

2. Computational Details 

The simulation was performed on 100 16^ x 32 
lattices at /3 = 5.6 with two flavours of dynamical 
staggered fermions, generated by the HEMCGC 
collaboration [3]. The light quark propagators 
were generated using the Clover action with tad- 
pole improvement, at three masses of light quark. 
For a comparison between using Clover and Wil- 
son light quarks, two of the clover kappa values, 
K = 0.1385 and 0.1401 were matched to the Wil- 
son kappa values 0.1585 and 0.1600 respectively 
using pion masses. Details of the Clover light 
spectroscopy results can be found in [4]. The 
heavy quark propagators were computed over a 
range of values of Mq between 0.8 and 10.0; the 
propagator in the static limit was also calculated. 
We will use a nommal value of = 2.0 GeV for 
this ensemble to present prehmmary quantitative 
results. 

We extracted the PS and V matrix element 
amplitudes in the standard way. We performed 
a simultaneous one exponential flt to SL and SS 
correlators, where the smearing function is the 
g.s. hydrogenic radial wavefunction with a radius 
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Figure 1. mgff for SL is offset by 0.2. 



of 3.0. Further details of our method and analysis 
can be found in [1]. The effective masses of the 
SL and SS pseudoscalar correlators for Mq = 
10.0 and k; = 0.1385 are shown in flgure 1; the 
mass obtained from the flt to the propagators is 
indicated. Note, even for a meson mass > 3Mb 
the signal does not fall into noise until t ~ 20, 
compared to t ~ 12 for the static case. 

The data at different Mq and k; are correlated 
and we found the optimal fltting range to be 7—20 
for PS and V correlators for k; = 0.1385 and 
0.1393 and 9-20 for ki = 0.1401, for allM^. For 
the static case we used 8 — 12. We computed the 
corrections to the currents separately by fltting 
to the ratio of the SL correlators with O2 and Oi 
at the sink, and in addition for the V meson the 
ratio of O3 and Oi . The perturbative calculation 
of Za and Zy is in progress, and has not yet been 
completed to include O3; we present the results 
for the tree-level operators. 



3. Results 



The PS decay amplitude with the tree level 



current correction included, Z^/vMpg , is plot- 
ted against the inverse meson mass (M2) in flg- 
ure 2 for K = 0.1385. The result in the static limit 
is also shown. We performed a correlated flt to 
the data using the functional form 
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Figure 2. 



range Q 


Co 


Ci 


C2 


C 1-8 0.2 


0.261(4) 


-0.25(1) 


0.11(1) 


W 1-7 0.7 


0.242(4) 


-0.25(1) 


0.13(1) 



Table 1 
Ki = 0.1385 



The fitting range was varied keeping the starting 
point fixed at the heaviest data point, Mq = 10.0. 
Beginning with a linear function, the fitting pro- 
cedure was repeated adding quadratic and cubic 
terms. In order to have confidence in the value 
for C'l a quadratic fit is required. 

We found it was possible to perform a linear fit 
to Za/VmIT for M^>3.0. A quadratic func- 
tion is necessary to fit to the rest of the data, 
ie there are significant 0(1/M^) contributions to 

/ corr . 

ZAjvMpi^ in the region of the B meson mass. 
Table 1 gives the values of C'i obtained from the 
best quadratic fit. The fit is shown in figure 2. 
The extrapolation to the static limit is consistent 
with the static result, as expected; the static case 
is the infinite mass limit of NRQCD. 

We compared these results with those obtained 
using Wilson light quarks at k; = 0.1585, also 
shown in figure 2. The coefficients obtained 

1 COTT 

from the best quadratic fit to Z^/vMps ^re 
given in table 1. There is a significant increase 

I COTT 

in Z^/vMpg |m°=oo of 10% using Clover light 
quarks compared to Wilson, and a similar de- 
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Figure 3. 



crease in Ci/Cq. 

Chirally extrapolating the Clover-NRQCD 
data and using the nominal value of a~^, we find 
Jb ~ 190(3) MeV and JbJJb ~ 1.18(3) (errors 
are statistical). This can be very roughly com- 
pared with js ~ 170 MeV and jsJjB ~ 1-3 for 
Wilson light quarks, using the same . 

A striking feature of the results is the large 
slope obtained; very roughly Ci/Co ~ —la~^ 
or ~ — 2 GeV, much larger than the naive ex- 
pectation of 200 — 500 MeV, and this leads to a 
large correction to the static limit at Mb of ~ 30- 
40%. In addition, previous lattice calculations us- 
ing Wilson or Clover fermions for the heavy quark 
around charm have found a significantly smaller 
slope, Ci/Co ~ — IGeV. To examine this further 
we obtained the individual contributions to the 
PS and V matrix element amplitudes from each 
of the 0(1/M) interactions. 

From eqn. (9) the linear part of the slope of 
{f\/Mps + 3/V^y)/4 arises purely from the ki- 
netic energy of the heavy quark. Figure 3 presents 
the results for ki = k^; renormalisation factors are 
omitted and this introduces a 0(as) ~ 20% error 
in the value of Ci/Co extracted. The PS matrix 
element amplitude with and without the current 
correction, which include contributions from the 
other 0(1/M) terms, are also shown for compar- 
ison. It is clear from the large mass region where 
the slope is linear that Gkin is much greater than 
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range Q Co Cj C2 



fVWp"' 1-6 1.0 0.227(9) -0.32(5) 0.23(7) 

SA 1-7 0.6 0.229(6) -0.29(2) 0.21(3) 

R 1-8 1.0 0.98(1) 0.33(10) -0.20(10) 

Rcorr ;^_g ;^ q 1.00(2) -0.18(9) -0.19(9) 

Table 2 
— ^^c 



the contributions from Guyp or Gcorr, and is the 
source of the large slope. The parameters ob- 
tained from the best fit to the data are shown in 
table (2) (SA). We find Gkin = C\/Co ~ -LSa"! 

or 2.6 GeV. 

To obtain the size of Guyp and Gcorr 
we consider the ratio of decay constants. 
fv^Mps/f^/Mv is presented in figure (4) both 
with and without the current correction. The 
data is consistent with 1 in both cases in the 
Mq = cxD limit, as expected. From a quadratic 
fit to all the data points we find Cp ~ — 0.2a~^ or 

0.4 GeV in table (2)). From eqn. (10) 

by omitting the current correction the linear co- 
efficient is determined purely by Ghyp- We find 
Ghyp ~ 0.03a-i or ~ 60 MeV (R in table (2)), 
and thus Gcorr ~ —O.Sa~^ or ~ —1.6 GeV. Since 
we calculate the correction to the current sepa- 
rately, Gcorr can be more accurately extracted 
from fitting to {02\PS) / {Oi\PS) directly. We 
find consistent results at this level of comparison, 

Gcorr ^ — 0.6(3 . 

A comparison can be made with the predic- 
tions of QCD sum rules. Using HQET Neubert [2] 

finds cp 2.9(5) GeV, cp„ ~ -0.9(l)GeV, and 

Gkin = —2.3(4) GeV, in good agreement with 
our results. As noted in section 1 Ghyp and Gcorr 
differ. In fact in HQET Gcorr = —A, and in con- 
tinuum NRQCD this becomes Gcorr = —Esim- 
Modulo lattice renormalisation factors we find 
agreement with this, Egim ~ 0.56. 

4. Conclusions 

We presented an analysis of the mass depen- 
dence of HL decay constants using Clover light 
quarks and NRQCD heavy quarks to 0(1/M). We 
find consistency between the extrapolation of the 
NRQCD data and the static case, and significant 
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Figure 4. Errors in the Mq = cxd extrapolation 
are indicated. 



0(1/M^) corrections to f\fM around Mp- The 
slope with I/M2 is significantly larger than pre- 
vious lattice calculations; this may be due to the 
use of too small a range of masses around the D 
meson when using Wilson or Clover heavy quarks. 
In addition, we obtained the three 0(1/M) contri- 
butions to the linear component of the slope sep- 
arately and found Gkin dominates. Good agree- 
ment is found between our results for Gkin and 
physical combinations of Ghyp and Gcorr and the 
predictions of HQET with QCD sum rules. The 
results indicate that the kinetic energy of the 
heavy quark gives rise to nonperturbative con- 
tributions to f\/M much greater than O(Aqcd), 
contrary to naive expectations. 

We found a significant increase in f\/~M com- 
pared to using Wilson light quarks; it will be in- 
teresting to see if improved scaling behaviour is 
found using a lower [3 value. 
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